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I. INTRODUCTION 

The topic of interacting ultracold fermions has wit- 
nessed some spectacular experimental developments 
since the first formation of a molecular Bose Einstein 
Condensate [T]. The Bose Einstein Condensate (BEC)- 
Bardeen Cooper Schrieffer (BCS) crossover has been the 
subject of numerous experimental and theoretical works 
[2] . Recently, the authors of Ref . [3] have found an accu- 
rate and elegant way of measuring the equation of state of 
an ultracold gas. In particular, they were able to measure 
the equation of state of a mixture of ^Li atoms in two 
internal states f and J,. In this experiment, the pressure 
P was measured as a function of the chemical potential fj, 
and the temperature T. The so called unitary limit (see 
below), where the interactions are known to be strong, 
was studied, starting from high temperature down to low 
temperature (close to the superfluid transition). In the 
high temperature regime, the data was compared to a 
virial expansion 



P{fi,T) = 2T{AT)-^z + b2Z^ + b. 



'3Z 



where 



solved p!5HT7] . we can calculate 63. 

The paper is organized as follows. In section |llj the 
general formalism is introduced. The technical key point 
is to work in (imaginary) time, rather than with frequen- 
cies. In section |III[ we apply the method to the deter- 
mination of hi and recover in a rather efficient way the 
Beth-Uhlenbeck result j[5j. Section IV is the main part 
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(1) 



of the paper, and it explains how to determine diagram- 
matically 63. Several technical details are given in the 
Appendices. The main analytical results are given in 
Eqs.([37l),(39),( |33| an d (43). Numerical results are shown 
in Figj?] and Eq.(j44[) . Finally, we conclude in section |v] 



II. GENERAL FORMALISM 

We describe now the method we develop to get an 
expansion in power of the fugacity. We first consider 
the Green's function for free fermions in imaginary time 
(for a general reference on diagrammatic techniques, see 
Ref. 18]). We have 



G°(p,r) 



{-e(r)+n^.(£p-/i)} (2) 



is the thermal wavelength (we take 



h = \ and the Boltzmann constant fc^ = 1). z = e^^ 
is the fugacity and /3 = T"^. In this way, the authors 
of Ref. [3^ were able to extract the numerical values of 63 
and 64. 

The theoretical calculation of virial coefficients started 
in the late 30 's with the work of Beth and Uhlenbeck 
[5] on the calculation of 62 (see also [3] and [S])- The 
problem of computing 63 was considered later in Ref. [7]. 
In the context of ultracold atoms, 63 for ultracold bosons 
was studied in Ref.jS] and later for ultracold fermions 
in Refs.[5HT^. It is the goal of this paper to present a 
diagrammatic method to calculate 63. Our results are in 
full agreement with Refs. fTU', and therefore disagree 
with Ref. [9] . The idea of the method is that one needs to 
include the 3-body problem in the manj/-body problem if 
one wants to calculate 63. In some sense, this is similar to 
the zero temperature low density expansion of Refs.[13l 
fT4] . In this work too, we will need to use the T-matrix 
of the 3-problem. Since the 3-body problem is basically 



with np{x) — (e'^^ -I- 1)"-^ and Ep — \i^/{2m) is the ki- 
netic energy of a fermion of mass m. Q{x) is the Heavi- 
side function. In the high temperature limit we consider, 
the fugacity z = e^^ is smaller than one , and we can 
expand the Fermi-Dirac distribution according to 



~/3e(p) 



1 + ze^l^^(p) 

\ ^ 1 \n-l. 



n>l 



z"(-l)"-^e- 



a/3£(p) 



(3) 



From Eqj2] and [3j we find an expansion in powers of the 
fugacity for the free fermions Green's function 



G°(p,r) = 



^G(°^")(p,r)z" 



ii>0 



(4) 



We have defined 

G("-°)(p,r) 
G(o^")(p,r) 



-e(r)e~^p^ (5) 
(-l)"-ie~"p^e-"'^"'', n > 1 (6) 
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Therefore G^"'") is retarded, while G^^'"'^ for n > 1, is 
not retarded. Notice that G'-^'"-' does not depend on the 
chemical potential /i. Diagrammatically, since G*^"'"-* is 
a retarded function, we represent it as a line with an 
arrow going from left to right if increasing time goes to 
the right (this is the Green's function of a particle in 
vacuum). On the other hand, G^^'"\ for n > 1, is not 
retarded, and we represent it as a n-times slashed line, 
which can be oriented from left to right or vice versa. In 
order to calculate the virial coefficients, we expand the 
density n{p,T) (per spin species for spin 1/2 fermions) 
in power of the fugacity z — e^^. 



n{fi,T) 



{ATr^z + 2b2Z^ +3b3Z^ 



nbr, 



+ ■ 



nip,T) = ^G(k,r = 0_) 



where we denote ~ J ,2^3 in 3 dimensions. 



n 



= ^G(°'i)(k,T = 0_) 
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(7) 



time 



FIG. 1: The lowest order diagram. 



III. CALCULATION OF 6, 



In order to calculate 62 F TOl: we must consider all the di 



which comes from Eq.Q and the Gibbs-Duhcm identity 
n = {dP/dfi)rp /2. The density per spin is related to the 
full Green's function G(k, t) via the equation 



(8) 



~ J (27r)3 

It is convenient to define the densities 6n^P^ such that 



(9) 



The principle of the method is the following : in order 
to calculate Sn'-P\ we need to find all the diagrams for 
G-t-(k, 0_) of order z^. Those are the diagrams with one 
G(O.P), or one G^^^P'^^ and one G^°-^\ or one G^^-p-^) 
and one G^^'^), or one G^^-p-^) and two G(°'1)'s etc... Of 
course, the higher the value of p, the more cases we have 
to consider, and the calculation becomes more and more 
difficult. Nevertheless, we have been able to use this 
method for the second and third viriel coefficients (62 
and 63). We notice immcdiatly that in order to calculate 
n(/Lt, r), hence G(p, 0_), all the e'^'^ terms in prefactors of 
the G^'s (c/. Eqjl]) cancel in the expressions. Therefore 
we do not consider them in the following. The chemical 
potential /i therefore appears only in the fugacity. 

At lowest order, we have simply G « G'^°'°^ + ^G(°'1). 
Qio,o) jg retarded and therefore does not contribute to 
G(k, 0_). Hence we find 



(10) 



grams with one G(°'2) (Fig 2(a) ) or 2 G^^'^^'s (Fig 2(b) I 



kt 

O 





time 



(a) 




time 



FIG. 2: The two diagrams contributing to the virial coefficient 
62. (a) is the diagrammatic expansion of the Fermi-Dirac 
distribution, (b) gives the Beth-Uhlenbeck contribution. 



The two diagrams contributing to ^n*^^-* are shown in 
Fig(2] The diagram of Fig |2(a)| gives the contribution 



= ^G(°'2)(k,r = 0_) 



(12) 



which is shown in the diagram of FigjT] The calculation 
is easily done using Eq|6] and we find 



,(1) - 



= (At)- 



(11) 



This is the ideal classical gas result. 

In the following, we focus on the problem of distinguis- 
able fermions of 2 species f (respectively I) of mass 
(respectively m^) typical of the BEG-BCS crossover (for 
a review, see [2]). Specifically, we consider equal masses 
{m^ = = m) and equal densities (unpolarized case) 
or equal chemical potential fx-f = /x^. 



One easily finds using Eqs|6] 



(13) 



This is of course just the free fermions contribution. 

We consider now the contribution of the diagram of 
Figj2jb). It contains explicitely the effect of interactions, 
through the 2-body T matrix. We show now that it gives 
the Beth-Uhlenbeck contribution [Sj.The imaginary time 
integration is in the time domain {0 < ti < T2 < /?}, 
corresponding to the domain for time differences {ti > 



3 



0, t2 > 0, /3 - {ti + t2) > 0}. We therefore find 



k,P • 



T2(P,t2)e- 



-(^-t2)(£tk+S4.P-k) 



(14) 



where e^,^k — k^/(2m^_^). The integration on ti is easily 
performed and gives a factor (/3 — t2)0(/3 — ^2)- The 
remaining integral on ^2 (for given P and k) is of the 
form 



dt2h{fi-t2)f2{t2) 



(15) 



with /i and /2 two retarded functions [i.e. which are 
equal to zero for a non positive argument). The convo- 
lution theorem for Laplace transforms gives 



/ = 



ds 



'0s 



Fiis)F2{s) 



(16) 



M — + = 2 TO 

the total mass of the 2-body problem, and (in 3D) 
, , 2n 1 

his) 



rrij. (a 1 — ^—2mrs) 



(17) 
(18) 



(19) 



This expression for t2{s) is easily found by writing down 
the integral (ladder) equation in time and taking the 
Laplace transform. This is just the usual result for a 
dimer propagator invacuum, with energy s. a is the scat- 
tering length and 



(m^ + m;) 



(20) 



the 2-body reduced mass. Next we go to the center of 
mass frame and make the change of variable k = k' + 
(771^ /M)P. By the change of variable s' = s- p2/(2M), 
the center of mass momentum P decouples. The inte- 
gration on P gives a prefactor (M/to^)'^/^(A^t) Here 
T is the thermal wavelength for f particles. The inte- 
gration on k' is also easily performed using (in 31?) 



E 

k' 



1 



1 



2tt \J—2mrS' 



(21) 



Finaly, we get the result 



M\ 



3/2 



ds' 
2jx~i 



t2is') 



to; 



■\/— 2toj.s' 27r 
(22) 



We show now that we recover the Beth-Uhlenbeck re- 
sult [S] from Eq. ( 22 1 . The integral on the variable 
s' is transformed by deforming the integration contour 
along the real axis. There comes two kinds of contribu- 
tions. The first contribution comes from the molecular 
pole of t2{s) (if there is one) at the molecular energy 
— = — l/(2TOra^). The second contribution comes 
from the branch cut along the positive part of the real 
axis (physically the continuum of scattering states). We 
finally get for the contribution to the density 



{rimol + nscatt.i 



C-y is the Bromwich contour of the problem. It is a with 
straight line in the s complex plane parallel to the imag- 
inary axis with a given real party 7. 7 is such that the 
integrand is analytic to the the left of the contour. We 
have also defined Fi{s) — J^°° dre''^/i(T) (and similarly 
for ^2(5)) for a complex variable s. The usual Laplace 
transform is defined with a minus sign in the exponential, 
but the usual expressions for propagators are recovered 
with this definition. 

In our problem, we have /i(t) = T'2(Pj''') and /2(t) = 
^e-^(^Tk-(-e4P-k), We easily find ^2(5) = l/(s - (etk + 
£;p-k))', and Fi{s) = <2(s - PV(2M)), with 



W = e^l^'^ifa > 
= 0, if a < 

+ CX3 



'^scatt. 



dxe-^"" 



y/2mrX (a 2 + 2mrx) 



(23) 

(24) 
(25) 

(26) 



We can recover the Beth-Uhlcnbeck [S] result by making 
the change of variable x — / {2mr) 



1 



+00 



rfpe-^^ ^ (27) 
(1 + p^a^) 



For equal masses, (to-|- = m^) we have for the second 
order virial coefficient 62 



&2 - 2 = ^2 {n,nol + riscatt.) 



(28) 



where Umoi is given by Eqs.([24|),([25|) and riscatt. by 
Eq.( |27[). I n the unitary limit, if a — > +00 we find from 
Eqs( |25|27l ) fijnoi — 1 and riscatt. = ^1/2 thus (for equal 
masses = to-|.) ^n^^'^'^Ai^ = \/2 (corresponding to 
62 = 1/^2 20J). If a -> ~oo, firaoi = and riscatt. = 1/2 
and we find the same result. 



IV. CALCULATION OF 63 

We now come to the calculation of the third virial coef- 
ficient &3 , which is the main result of this work. We focus 
on the 3-D situation but the method works in any dimen- 
sion. We proceed following the general method presented 
in section [IT] and developped in section |III[ We want to 
gather all the diagrams of order z^. There are three sorts 
of diagrams (the partition number of 3) : 1) with one 
G("'3)^ 2) with one G^^'^) and one G(°'1), 3) with three 
(5(0,1) 'g^ We will denote the corresponding corrections 



4 



to the density Sn^^'^\ Jn^'^'^' " ^ and 5n'^^'^'"''> (there are 
several contributions in the last two cases). 

The diagram 1) is nothing but the free fermion result 
and it amounts to expand the Fermi-Dirac distribution 
in power of z. This is a standard result, and we find 
for the density the correction for free fermions (in d = 3 
dimensions) 



fc(3'i)A3, = 3- 



■3/2 



(29) 



similar to the diagram of Figpjb). Indeed, by noticing 
thatG(0'2)(k,-(ti+t2)) = (-e-/5-'')G("-i)(k,-(ti+t2)) 



we just have to multiply the time integral in Eq.(14) by 
_g-/3etk in order to get the contribution of diagram of 
Figjsf^a). We can therefore skip directly to Eq.([l6| in 
order to evaluate the time integrals. Then we make the 
following change of variables (the Jacobian is unity) 



The diagrams with one G^°'^^ and one G'-*^'^-' are shown 
in Figj3] The diagrams of Fig |3(a) and |3(b)| are very 



kt 



P-kj 



T2 




kt 



kt 



-time 



P-ki 



T, 




kt 



-time 



k' 
P' 



4m 
2 

(k + P) 



(30) 
(31) 
(32) 



(a) 



(b) 



FIG. 3: The two diagrams contributing to 63 with one 
and one G(°'^' (see text). 



The integral on P' decouples and can be done analyt- 
ically. The angular integration on k' can also be per- 
formed, and we get 



3^/3 



+ CX3 



dk' k''- 
27r2 



ds' 



his') 

fs'- 



(33) 



The diagram of Fig |3(b)| can be calculated along the 
same line. We just have to exchange the role of f and I 
particles. As a result we find 

At unitarity (a^^ = 0), the integrals can be done ana- 
lytically. 



_ / 8 2 
ViVi " 3^^ 



B. 5n(^'3) 



'0.300994 (35) 



diagram also always contains an outgoing k, f line from 
the imaginary time r = 0. This can be a G(°'i) (slashed) 
line (first case) or a retarded G^^'^^ line {second case). 



1. 5n(3.3.i) 

In the first case, the only possible diagram is shown in 
Figj4]and is similar to diagrams of Figs |2(b)| and Figjs] 
If we follow the same line of reasoning than for the cal- 
culation of these diagrams, we easily find 



(36) 



We want to find all the diagrams with three G'-'^'^-' 
propagators. A given diagram always contains an incom- 
ing k, t line to the imaginary time t = 0. The initial time 
of this fermionic line must be in the interval [0, /?] . There- 
fore, this line can not correspond to a retarded propaga 



2. 5n(3'3'2) 



The second case is the most interesting. A general 



tor G^"''^', and it must be a Gj"'^"* (slashed) line. A given diagram is drawn in Figjsj It contains two G^'^'^^s and 



5 



P-ki 




T2 








time 



FIG. 4: The diagram contributing to 5n'^'^'^' 




FIG. 5; The general structure of a diagram contributing to 

5„(3,3,2) 



one G 



(0,1) 



(Fig 



(Fig|6(e)]and[6lT)^ 



6(a) to 6(d) I or two g|°'^'''s and one G 



^(0,1) 



We consider first the case with two g1"'^'''s and one 



G 



(0,1) 



The incoming \, (from the left) can first inter- 



act with the incoming (from the left) G^°'^'' (diagrams of 
Fig' ' ' ' 



6(a; 



and 



6(b) ) or with the incoming g'^'^^ (diagrams 



of Fig |6(c) and |6(d) ) . In each case, the two diagrams are 



found by choosing which line G^^''^^ the outgoing (to the 

right) 4, can last interact with. In all these diagrams, Tg 
is the 3-body T- matrix, extensively used in [131 El [H! ■ 
It physically describes the scattering of a f particle with 
a dimer (diagrammatically represented by T2). 



Gf^)'s 



The case with two 
treated in a similar way 
line interact with the same G 



and one Gj ' can be 
in Fig |6(e) the external f 

''"^'^^ line at the entrance 
4. 



(on th e left ) and at the exit (on the right) of Tg ; whereas 

's have been exchanged at the 



6(f) 



the two g|°'^^ 



in Fig, 
exit of Tg^. 

We first consider the calculation of 5nS^'^''^''^\ This 
is explained in Appendix]^ The difficulty is to perform 
the time integrals. It is however possible to transform a 
multidimensional time integral into a simple integration 
on a (complex) energy. Mathematically, this comes from 
the fact that the time integrals are just convolutions, and 
the convolution theorem for Laplace transforms enables 
to write it as an inverse Laplace transforms of products 
of functions of the energy, as we did for the calculation 
of 62- We therefore give the results 



^^(3,3,2,a)^3^ = 



2tt ./c 27r i 



+00 ^^^2 



dpp 



t2{s' 



27r2 



-2 mr{s' 



2 rriAt.J 



-) \l>0 



(37) 



C is a Bromwhich contour such that the integrand is an- is the atomf-dimer reduced mass, 
alytical to the right of C. 

r : \ o The calculation of Sn'-'^-'^'^'''^ is similar, and we just give 

m+(m-t- -t- mi ) 2 j j t> 

"T-At D = -ji ^ = (38) the result 

(2m^ + m^) 3 



^^(3,3,2,6)^^3^ = 3^/2 V / ^e-f^^' f 

,-^r, Jc i Jo 



i>Q ■ 



PP 



+°° dp 
2^ 



P%{s' 



2 m 



At,D Jo 



dp' 

2^ 



p"h{s' 



P' 



-) 



(39) 



r 



with Qi{z) 



dQjS^ _ 1 jvL^ ^y_ii(^ i^Q^ jg ^}jg Leg- endre function of the second kind). 



6 




(e) (f) 
FIG. 6: The 6 diagrams contributing to (5n'^'"^'^^ (see text). 



The numerical result is shown in Figj?] In order to reveal 
the effect of interactions, we have plotted the difference 
63 — 3^^/^ as a function of At /a. In practice, in order 
to integrate safely on s', we deform the contour C into 
two semi straight lines, symetric with respect to the real 
axis, and making an angle a G]0,7r/2[ with it. In this 
way, the exponential goes to zero at the extremities of the 
contour. We checked for different values of a, that the 
result is independant of a. We discretized the values of s' , 
p and p' on a grid in order to determine ^ numerically 



from the linear integral equation ( B7 1 , which transforms 



into the problem of solving linear algebraic equations. 
Finally, for a given value of we sum on s',p and p' 
in order to get i5n(^'^'^'°) and Jn^^'^-^'''). The curve of 



CO 

I 




-1.5 



-3 -2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 

At /a 



The calculation of diagrams of Fig 6(c)| and |6(d)| is 
done in the same way and we find that 



5„(3,3,2,c) ^ ^„(3,3,2,6) 



(40) 
(41) 



Moreover, it is also easy to verify that Jn'^''^'^''^-' (re- 
spectively Sn^^'^'^'-f^) is obtained from the expression for 
j^(3,3,2,a) (respectively Sn^^'^'-^'''^) by exchanging the role 
of t and J, particles (this amounts to exchange and 
nil analytical expressions). 



C. Calculation of 63 in the case rrif — mi 

In general, 63 is obtained from Sn^^^^ by summing all the 
contributions ^^(3,1) ,5yj(3,2,aandb)^ j^(3,3,i)^ Q^Yi^ the six 
^j^(3,3,2)5g^ However in the case of equal masses, we have 

J„(3,3,l) ^ J„(3,2,a) ^ g^(3,2.b) ^ ^„(3,3,2,a) ^ J„(3,3,2,d) ^ 

j,^(3,3,2,e) ^ ^Yiis way we get 



FIG. 7: The third virial coefficient 63 
the dimensionless parameter At /a. 



- 3 ^^■^ as a function of 



FigjTjis in excellent agreement with the results of Ref . [TU] . 
At unitarity, we find for the 63 coefficient (we sum the 
harmonics in Eqs.(37),(39 1 up to Z = 7 for the four digits 
we show) 



-5/2 



-0.3551 



(44) 



in excellent agreement with the results of Refs[TOl [12]. 
This is therefore an independent check for the value of 
&3, since our method is very different from the method of 

[101 [ni. 



V. CONCLUSION 

We have calculated the third order virial coefficient 
&3 in the problem of 2-species fermions interacting via 



3/2 ^_ 3 1^^^(3,2,0) _j_ ^^(3, 3, 2, a) ^^(3,3,2,6)\ ^3^a short range interaction. We have developped a dia- 
^ ' grammatic method which leads to explicit analytic ex- 

(42)pressions. This approach might be extended to the cal- 
culation of the 4th order virial coefficient 64, which was 
recently calculated 1121. 

Acknowledgments: We thank R. Combescot, F. 
Chevy, C. Mora, N. Navon, N. Regnault and F. Werner 
for stimulating discussions. 



And therefore for 63 = (5n(^)Af,/3 

63 = 3-5/2 ^ |^J„(3,2,a) ^ ^„(3,3,2,a) ^ ^^(3,3,2,6)^ ^3^ 

(43) 
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Appendix A: Calculation of (5n(^'^'^'"' 



In order to calculate the diagram of Fig 6(a) we first 



look at the "Born" approximation of tJ. The diagram 
is shown in FigjSj In this figure, we have indicated the 
4 relevant time differences ii,...4 as well as a convenient 
way to choose the momentum variables. Physically P is 
the total momentum of the 3-body problem (2 f, 1 l). 
p is the momentum of the incoming f, in the centre of 
mass reference frame, scattering with a dimer of opposite 
momentum in the centre of mass reference frame. We 
have defined 



a 



{2m^ + m^) 3 
We give expressions valid for any masses and m-^. 



(Al) 



The time domain is terms of time differences is {^i > 

0, • • • ,t4 > 0, /3 - (^1 H 1-^4) > 0}. First, we notice 

that the integrand is independent of <i. Integrating on ti 
gives a factor (/3 - (ia + ^3 + ^4)) B (/3 - (^2 + ^3 + ^4))- 
After integration on ti, we can then write this diagram 
as (we do not sum on momenta, and denote it /) 




I = - J ^t=2dU (T2(-p + (l-a)P,t2)e-^'^(P+"P)*^e(t2)) (e-(2ei'(P+"P)+^^(-2p+(i-2a)P))t30(i3)^ 

X (t2 i-p' + (1 - «)p, u) e-^^(p'+"P)*^e(t4)) 

X ((/? - (<2 + i3 + h)) e (/3 - (<2 + ^3 + h)) e-(/3-(*2+t3+t4))(et(p+oP)+et(p'+aP)+e,(-p-p' + (l-2a)P))^(^2) 



This apparently complicated expression is actually a We then easily find 
convolution integral of the type studied in Appendix [C| 



ds 



-fis 



t2 S 



2(2rn^ + m;) 2mA-t,D 



P2 

2(2m-|-+mj^ 



0-i 



(<' 



p2 



2{2m-x+iyn) I 2 V m-f ^ mj, ■ 



(A3) 
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If we look at the next order Feynman diagrams for T3 , we 
see by inspection that we still find a convolution integral, 
and we just have to replace the last line of the above 

equation by 4(P>P; " 2(2m^+m^) ))- ^^^^ ^^a^' ^he 
multimensional time integrals are replaced by a single 
integral on the variable s. 

We now come to the summation on wave vectors. For 
the summation on the total momentum P, we see that, 
similarly to what was found in the calculation of 62, P 

enters only through the combination s ~ „,„ ^ . This 

enables to decouple the summation on P (by the change 
of variable s' = s — , ^ ). After summation on P, 

3/2 



2(2m-|-+mj^) > 

3 / 2m-f+7ni ' 



we get a factor {JS-i^t) 

Finally, for it is also possible to inte- 

grate on p', since p' enters only the expression through 
7^ ^ — = — r^2^- The integration is done 

via the change of variable p" = p' + (m^)/(77i-|- + m^)p, 
and the use of Eq.(21|. In this way we get a factor 



(m,)V(27r)/V-2TO.(s'-pV(2TOAt,D)- 

The last step is to express i3(p,p; s') in terms of the 
Legendre Polynomials components t3;'s. This is done 
thanks to 



4(P'P; 



(A4) 



In order to be as general as possible, we also give her e 
formulas valid for any masses and m^. In Eq.(Bl|, 
it is possible to perform the integration on the variable 
Sq by deforming the integration contour in the half plane 
Re{sq) > 7. The only singularity in this domain is the 
pole Sq = e-fq, and we end up, as usual, with an integral 
equation where Tg in the integral term is evaluated "on 
the shell" . The next step is to make a change of variables 
and function according to 



Pi 

P2 



Pi 
P2 

Sp 



2 + 

"^T p 

2 m-f + 

P^ 

2(2 m^ + m^) 



(B2) 
(B3) 

(B4) 



4(Pl.P2;Sp) = rl({pi,£^(pi)},{p2,£t(P2)};{P,Sp}) 

(B5) 



which physically represents the forward scattering of an 
atom on a dimer (we used P;(l) = 1). In this way, we 
recover Eq. ( 37 1 . 



Appendix B: Integral equation for the 3-body 
problem in the center of mass frame of reference. 



We first recall the integral equation for Tl{pi,p2;P) 
(for complex frequencies) following the approach of [T7] 

- Gt{q)G^{P -pi~ q)T2{P - q)Tliq,p2;P) 
q 

(Bl) 



This amounts to go to the center of mass reference frame 
and to evaluate T3 "on the shell" . p'l is the momentum of 
a t incoming atom in the center of mass reference frame 
(total impulsion P in "laboratory" frame), s'p can be 
seen as the total "energy" minus the center of mass ki- 
netic energy, i. e the energy in the center of mass reference 
frame. In this way, we get 



4(pi-P2;s'p) 



{p'i)^+ip'2V I p'rP2 

2m,, mj^ 



All 



(27r)3 



2mr ni^ 



4(q':P2;sp) 



(B6) 



We do not suppose here that the two species have equal 
is the 2-body reduced mass, and 



masses, 



(m^H-mj.) 

WAt,r> — "(2m^+m^j"' ^® ^^"^ atom^-dimer reduced mass. 



This integral equation can be projected on the Legen- 
dre polynomials Pi. In this way, the integrals equation 
for different I's decouple, and we find 
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/ \2 I /„/\2 



><tUq\P2-,s'p) (B7) 



where Qi{z) — ^ J^^du^^^Pi{u) is the Legendre func- 
tion of the second kind 121, . 



Appendix C: Convolution integrals 

We want to compute the following 4-diniensional inte- 
gral 



/(/?) = / Tlt2dU0W-it2+h+U))g{P-{t2+h+U))f2{t2)f3ih)Mt4) (CI) 
Jti>0 

I 



Wc first integrate with respect to ^2, and use the convo- 
lution theorem for Laplace transforms 

f dt29{t~h)f2{t2) = C-^ {G{s2) F^is^)) (C2) 
Jo 

where G{s) = J^°° dte''g{t) and ^2(5) = 

dte"^ f2{t). We use a notation such that the 
Laplace transforms of a given function is denoted by 
the corresponding capital letter. The symbol C denotes 
the Laplace transform, and C^^ the inverse Laplace 
transform. The integration with respect to t2 therefore 
gives the contribution 

g2{P^{h + U)) = /:-MG(s)i^2(s))*=^_(,3+4C3) 

The integration on ^3 can be done similarly, and we get 

53(/3-i4) = C-^{G2{s)F3{s)\^p_,^ 

= {G{s) F2{s) F,{s)\.,^ (C4) 



The last integration on t4 gives the full integral 

m = C-'{G3{s)F,is)),^^ 

= (Gis) F2is) F,is) F4s)\^p (C5) 

In this way, we have transformed the 4-dimensional inte- 
gral Eq |C7| in the 1-dimensional following integral 



I{I3) = / ^e-f'^G{s)F,{s)F,{s)F,{s) (C6) 

with Cj a contour parallel to the imaginary axis with 
Re{s) — 7, such that G,F2,F3 and F4 are analytical for 
Re{s) < 7. 

This can be generalized easily to any n-dimensional 
integral of the form 



/ Ii1^2dU0{l3-{t2 + ---+t^))9{^-{t2 + ---+t^))f2{t2)---fn{tn) (C7) 

Jti>a 

I 

Following the same method as before, we get 

/„(/?)= / ^e-^^G{s)F,{s) ■■■ F^{s) (C8) 
Jc-, 27ri 
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